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ABSTRACT 

Surface brightness profiles for 23 M31 star clusters were measured using images from 
the Wide Field Planetary Camera 2 on the Hubble Space Telescope, and fit to two types 
of models to determine the clusters' structural properties. The clusters are primarily 
young (~ 10 8 yr) and massive (~ 10 4 ' 5 M Q ), with median half-light radius 7 pc and 
dissolution times of a few Gyr. The properties of the M31 clusters are comparable to 
those of clusters of similar age in the Magellanic Clouds. Simulated star clusters are 
used to derive a conversion from statistical measures of cluster size to half-light radius so 
that the extragalactic clusters can be compared to young massive clusters in the Milky 
Way. All three sets of star clusters fall approximately on the same age-size relation. The 
young M31 clusters are expected to dissolve within a few Gyr and will not survive to 
become old, globular clusters. However, they do appear to follow the same fundamental 
plane relations as old clusters; if confirmed with velocity dispersion measurements, this 
would be a strong indication that the star cluster fundamental plane reflects universal 
cluster formation conditions. 
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Introduction 



The spatial distribution of stars within a star cluster is an important indicator of the cluster's 
dynamical state, and the structural parameters (e.g. core, half-light, and tidal radii; central surface 
brightn ess, and concen tration) indicate on what timescales the cluster is 'bound' to dissolve. The 
work of Soitzer ( 198 71) showed that core collapse is an inevitable part of cluster dynamical evolu- 



tion. IDjorgovski fc King! (119861 ) were am ong the first to determine t he fraction of core-collapsed 
Milky Way globular clusters (GCs), while iDiorgovski Mevlanl (| 19941 ) examined a large sample of 
Milky Way clusters and defined the 'fundamental plane', showi ng that surface brightness pr ofiles 
of Galactic GCs were well-described by only a few parameters. iMeylan fe Djorgovskil (119871 ) sur- 
veyed GCs in the LMC and SMC for core collapse and found that only a handful of clusters were 
core-collapse candidates; they suggested that environmental or age effects were responsible for the 
difference with Milky Way globulars. 

A few s patially - resolv ed st udies of GCs beyond th e Magellanic Clouds were done with ground- 
based data. iRacinej (|199ll ) and iRacine Harris! (I1992l) used high-resolution i maging to distinguish 



M31 GC candidates from background galaxies, and lCohen &: Freeman! (|199ll ) determined the tidal 
radii of 30 M31 halo GCs, finding them to be similar to Milky Way GCs. However, detailed studies 
of the structures of M31 G Cs awaited the angular resol ution of the Hubb l e Spa ce Telescope. The 
first work on M31 GCs by iBendinelli et al.1 (119931) and iFusi Pecci et al.l (119941) was follow ed by 
numerous oth e rs inc l uding Rich et al. (1996), Grillmair et al. ( 19961 ). Holland et al. (1997), and 
Barmbv et al.l (|2002l . 120071 ). Clusters in Local Group galaxies are near the limit for resolution 



into individual stars by the Hubble Space Telescope (HST), although some structu ral information 



such as half-light radii can be recovered for clusters in more distant galaxies (e.g., lHasegan et al 



20051 ) . Conclusions of the studies of extragalactic globulars incl ude the dependence of cluste r size 
on galactocentric rad ius, first pointed out for the Milky Way by IDjorgovski &: Mevlanl (119941 ) and 
van den Berghl (|1994l ): a pos sible differenc e between sizes of clusters in different metallicity groups 
(for a detailed discussion see IJordanl 120041 ): and a recognition that globular clus ters in a variety of 
environments appear to lie on the same fundamental plane (jBarmbv et al.l 120071 ) . 



Structural studies of younger star clusters present more difficulties. Open clusters (OCs) in 
the Milky Way are generally much less massive than globular clusters. As viewed from our location 
in the Milky Way, they are embedded within the disk, so that the cluster is easily lost against the 
much more numerous field stars, and determining stellar membership in these less-concentrated 
objects is not straightforward. Comprehensive studie s of Milky Way open clusters are relatively 



recent: 



Kharchenko et al.l (|2005l ) and follow-up work (jSchilbach et al.l 120061 : iPiskunov et al. 



2007 



20081 ) measured a variety of radii (core, c orona, tidal) for s everal hundred clusters and found their 
masses to be in the range 50-1000 Mq. iBonatto Bical (|2005l ) analyzed in more detail a much 
smaller number of Milky Way open clusters, finding that the cluster size increased with both age 
and Galactocentric distance. These authors also found that their sample of clusters showed evidence 
for an 'open cluster fundamental plane.' 
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Milky Way open clusters are not the only known population of young star clusters, and possibly 
not even the best one to study. The Galactic OCs cover a limited range in age and mass and their 
census is suspected to be far from complete because of extinction in the Ga lactic plan e. The 



Magellanic Clouds (MCs) have many young star clusters, recen tly cataloged by iBica et al.1 (|2008l ). 
The brighter MC clusters were studied in a pioneering work bv lElson et al.l (|1987l ). These authors 
analyzed the radial profiles of 10 clusters and found t hem to be b etter-fit by 'power-law' profiles 
of the form I(R) oc [1 + (i^/ rn) 2 ] - ^ 7-1 ^ 2 than by the King J 1960) models conventionally used to 
fit globular cluster profiles. McLaughlin van der Marell |2005h re- analyzed a large set of MC 
cluster data and found the situtation to be somewhat more complex. Those authors argued that 
the extended envelopes characteristic of the power-law profiles are a generic feature of many young 
and old star clusters and that "the development of a physically motivated model accounting for 
this . . . could lend substantial new insight into questions of cluster formation and evolution." 

Outside the Milky Way, many galaxies are fo und to have 'young massive clusters' (YMCs; 



Holtzman et al 



(Brodie et al 



1992; IWhitmore Schweizerl Il995l ). These clu sters have ages up to a few Gyr 



1998 ) and masses comparable to globular clusters ( Larsen fc Richtlerl 1999 ). Studies 



of YMC structures show correlatio ns of power-law slope 7 with age (jLarsen 



rate (Weidner et al 



20041 ). core radius 



with age (Mackcy &: GilmoreJ [2003J) , and mass of the brightest cluster with galaxy star formation 



2004 ). As of yet there is no comprehensive study of star cluster structures 



over the full age and mass ranges seen in nearby galaxies. M31 is now recognized to also have 



a large population of young star clusters (|Fusi Pecci et al. 



2005 



Caldwell et al.ll2009h . although 



their relationship to both the YMCs and globular clusters is not well-understood. The purpose 
of this paper is to carry out an initial study of the structural properties of some young M31 
clusters. We analyze a sample of 23 clusters using data from the Wide Field Planetary Camera 
2 (WFPC2) onboard the Hubble Space Telescope; extensive analysis of 'artificial clusters' (see 
Appendix) i nforms our analysis proced ures. Throughout this work we assume a distance to M31 
of 783 kpc (jStanek Garnavichl Il998l ). for which 1" corresponds to 3.797 pc. All magnitudes 
are in the Vega sy stem, and cluster names use the convention of the Revised Bologna Catalog 
(|Galleti et alJboO^ I 2 ! see that work for cluster coordinates and other properties. 



2. Data and analysis methods 
2.1. Cluster sample 



The study of star clusters in M31 has a long history dating back to at least iHubbld (|1932l ). so 
any attempt to assemble a sample of young massive clusters necessarily draws on many previous 
works. While a number of studies of the glo b ular cluster system have noted the presence of possible 
young clusters in M31 (IBarmby et al.l 12000 : IWilliams Hodgdl200ll ). the first comprehensive list 



2 Online version at \protect 



lttp:/ /www. )o.astro.it /M31 
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of such objects was assembled bvlFusi Pecci et al.l (120051), who called t hem 'blue luminous compact 
clusters', or BLCCs. krienke & Hodeel (I2OO7I I2OO8I) andlHodee et al.l (|2009l ) searched for M31 'disk 
clusters' in archival HST imaging data, and I Caldwell et al.l (|2009l ) pr esented a comprehens ive list 
of nearly 150 young cluster candidates from a spectroscopic survey. Caldwell et al] ( 2009) noted 



that the handful of their young clusters with measured structural properties (from iBarmbv et al 



20071 ) covered a wide range in p arameter space. Th e HST resolved-star study of four 'massive and 
compact young star clusters' by lWilliams &: Hodgd (|200ll ) (program GO-8296) did not include an 
analysis of the objects' structural properties. 

The ma i n samp le of clusters studied here is described in detail by the companion papers by 



Perina et al. (2009 



32 



The present project began with an interest in confirming the results of 
Cohen et al. (120051 ) who used adaptive optics imaging to s how that some of the clusters proposed 
as young were in fact a sterisms (but see the contrary view of lCaldwell et al.ll2009l and the discussion 



111 



Perina et al 



2009 a|). HST progr am GO-10818 was aimed at imaging all of the 'class A' clusters 



proposed by (IFusi Pecci et al.ll2005l ) which did not already have HST imaging, a total of 21 objects. 
In the course of the program we found that two clusters in the candidate list were in fact the 
same ob ject (IPerina et al.1 l2009al ) . and the object NB67 was a star, so the program contains 19 
objects. IPerina et al.1 (|2009bl ) showed that 16 of the clusters are young, with ages < 1 Gyr, and five 



(B08 3, B222, B347, B374, and NB16) are in fact intermediate-aged or old (see also ICaldwell et al 



2009). We retain these five clusters in our sample but show them with different symbols in the 
analysis. We augmented the GO-10818 data with archival data on the four clusters studied by 
Williams fe Hodgd (|200ll ) to bring the total number of clusters to 23. HST archival data exists for 
additional clusters but in the interests of dealing with a mostly-homogeneous dataset we restricted 
the sample to only the GO-10818 and GO-8296 clust e rs. Th ree of the clusters in the latter dataset 
had structural parameters reported in lBarmby et al.1 (120021 ): here we re-analyze them in a manner 
consistent with the other clusters. Except for B083 a nd B347, all of the clusters are projected 
against the M31 disk (see Fig. 1 of lPerina et al.ll2009bl ). 



2.2. Data reduction and surface brightness profiles 

The GO-10818 program was originally intended to be carried out with the Advanced Camera 
for Surveys (ACS), but because that instrument failed, the images were obtained instead with the 
Wide-Field Planetary Camera 2 (WFPC2). All objects were observed with two 400-s dithered 
images in ea ch of 2 filters: F450 W and F814W (for further detail, and an example of the CMD 



analysis, see IPerina et al.ll2009al ). The GO-8296 program was also carried out with WFPC2 and 
involved two 800-s images in F439W and two 600-s images in F555W (as well as longer images in 
F336W which are not used here). The target clusters were on the PC chip in all cases, and only 
data from that chip is used in the present analysis. Table Q] summarizes the datasets together with 
other pertinent information about the clusters. 



The multiple images were combined with the STScI Multidrizzle software, using the 'recipes' 
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provided on the drizzle webpage. The pixel scale of the resulting images was 0.0455", or 0.172 pc at 
the M31 distance. While correcting for Charge Transfer Efficiency losses would be desirable, there 
is currently no prescription available for correcting surface photometry of extended objects so no 
correction has been made in the present analysis. Although M31 star clusters are relatively large 
(a few arcsec) compared to the HST optical point-spread function (PSF), convolving model profiles 
with the PSF prior to comparison with the data should improve the accuracy of measurements of 
the cluster cores. Model PSFs were generated for the relevant filters at the camera center using 
TinyTim. The clusters are small compared to the camera field-of-view, and PSF variation over the 
cluster extent is negligible. 

Transform ing instrument a l mag nitudes to calibrated surface brightness was done following the 



prescription in iBarmby et al.1 (|2007l ) . Image counts were first multiplied by the inverse square of 
the pixel scale to give counts C in units of s -1 arcsec -2 . These can be transformed to magni- 
tudes arcsec -2 through \i = Z — 2.51og(C), where Z is the instrument zeropoint. They can also be 
transformed to intensity I in L & pc -2 through / = 10 0A ^ Z '~ Z ^C . (Independent of the instrument 
used, Z' = (m — M)m31 + M® + 51og(/3) = 21.5715 + M® where j3 is the number of arcsec corre- 
sponding to 1 pc; (3 = 0.2644 at the assumed distance of M31.) The zeropoints used come from the 
respective instrument handbooks; the solar magnitudes are from calculations by C. WillmeJl. All 
are listed in Table [2] for reference. 

Studies of surface brightness profiles of Local Group star clusters are in a somewhat different 
regime from either Galactic clusters or clusters in more distant galaxies. Local Group star clusters 
are resolved into stars in their outer regions but not in their cores. They differ from galaxies 
with comparable angular sizes (< 10 arcsec for M31 and M33 clusters) in that the galaxies are 
composed of many more stars and have much smoother light distributions. To better understand 
the limitations of our analysis, we simulated artificial star clusters, measured their surface brightness 
profiles, and fit those profiles to models: these simulations are described in Appendix lAl 

Surface brightness profiles for the M31 clusters were measured by combining integrated pho- 
tometry with star number counts (the 'hybrid' procedure described in Appendix |A|) . In the inner 
regions of the clusters, surface brightness profiles were derived using the IRAF ellipse package 
to fit circular isophotes to the image data. The isophote centers were fixed at a single value for 
each cluster, with centers determined as the intensity-weighted centroid in a 75 by 75 pixel box. 
Star counts were derived only from stars within specified regions of the CMD, with the designated 
region varying by clu ster depending on th e age. The details of the star counts for the GO-10818 



clusters are given by iPerina et al.l (|2009bl ) ; for the GO-8296 clusters, star counts were computed 



from background-subtracted CMDs (Fig. 6 of lWilliams &; Hodgdl200ll ) with positional data kindly 



provided by B. Williams. The star counts were used for radii > 7 pc (40 pixels) from the cluster 
centers, and scaled to linear intensity units (L® pc -2 ) by matching the counts and photometry 
over the overlap region 5-10 pc. The same star counts were matched to integrated photometry 



J http: //www. ucolick. org/~cnaw/sun. html 
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profiles in both red and blue filters, but with different scaling factors; star count uncertainties were 
matched to the photometry uncertainties by scaling as for the intensity. No background subtraction 
was performed on the star counts. 



2.3. Profile-fitting methods 



There are a number of po ssible choices for st ar clus t er de ns ity profiles, i nclud ing King! ([1966, 
hereafter King), IWilsonl (119751. hereafte r Wilson), iKingj (|1962l ). lElson et al.l (|1987l . also known as 
'power-law' or 'EFF'), and ISersid (|1968l ). Unlike the other three types of model profile, the King 
and Wilson models have no analytic expressions for density or surface brightness as a function 
of projected radius; profiles are obtained by integrating phase-space distribution functions over 
all velocities and then along the line of sight, assuming spherical symmetry (for a review, see 
McLaughlin! 12003) . The King model is t he most commonly-used in studies of star clusters; how- 



ever, 



McLaughlin &: van der Marell (|2005l ) showed that, with data that extends to sufficiently large 



projected radii, many Local Group clusters are better-fit b y the more-extended W ilson models. 
Globulars in NGC 5128 are also better-fit by Wilson models (McLaughlin et al.ll2008i ). although an 
analysis using nearly identical techniques (IBarmby et al.l 120071 ) found that massive M31 glob ulars 
were b e tter-fit by K i ng mo dels. Taken t ogeth er , these recent analyses showed that fitting the King 
(|1962l ). lElson et al.1 (119871 ). and ISersid (119681 ) models did not add significant information beyond 
that provided by the King and Wilson models, so we consider only these two models in our analysis. 

The King and Wilson models are single-stellar-mass, isotropic models defined by phase-space 
distribution functions of stellar energy E: 



f(E) oc ^ 



exp[ 
exp[ 




-E/al 
-E/al 



1 + E/al 



E < 
E < 
E > 



(King) 

(Wilson) 

(both) 



(1) 



where o"o is the central velocity dispersion. The effect of the extra term in the Wilson model f(E) 
is to make clusters more spatially extended. Both sets of models are characterized by three param- 
eters: a dimensionless central potential Wo, which measures the degree of central concentration; 
a scale radius ro, which sets the physical scale; and a central intensity Io, which sets the overall 
normalization. For the King models, Wo has a one-to-one correspondence with the more-familiar 
concentration c = log(rt/ro), where rt is the tidal radius at which the density p{rt) = 0. Possibly 
contrary to intuitive expectations, for two profiles with the same scale radius, the profile with a 
larger value of c or Wq declines more slowly. 

Deriving the structural properties of the simulated clusters involve d fitting their pro jected sur- 
face d ensity pro files to models using the GRIDFIT program described bv lMcLaughlin &: van der Marel 
(|2005l . see also iMcLaughlin et al.ll2008l ). The program uses a grid of model density profiles, pre- 
computed for a range of values of Wo, then finds the scale radius ro and central surface brightess 
Io to minimize the weighted x 2 f° r each Wo) the best-fitting model is the one with the global % 2 
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minimum. The model profiles are convolved with the instrumental PSF before comparison to the 
data. Since no background subtraction was performed on the star counts, the background level was 
determined as one of the parameters of the model fitting. For a few clusters the fitting algorithm 
converged to unreasonably large or small values, and a fixed background corresponding to the 
lowest level reached by the star counts was subtracted before re-fitting; in general this procedure 
improved the reduced x 2 °f the fits. 



2.4. Profile- fitting: results 



Figure Q] shows the cluster surface brightness profiles together with the best-fitting models. 
The parameters of the models are give n i n Table El corrected for e xtinction using the values of 
E(B — V) given bv lPerina et al.l (|2009bl ) or lWilliams fc Hodgd (|200ll ). Conversion of filter-specific 
measurements to the F-band is done using the transformations described in the appropriate HST 
Instrument Handbooks; briefly, we compute the extinction-corrected color (V — x)o, where x is the 
observed-band magnit ude, as a fun c tion o f color in standard bands (e.g., (V — I)o)- Ground-based 
integrated colors from lGalleti et ak ( 2007 ) are used for the standard-band colors, to avoid iteration; 
uncertainties of 0.1 mag in (V — x)n are assumed and propag ated through the parameter estimates. 
As previously shown by lMcLaughlin fc van der Marell (120051 ). differences between Wilson and King 
model profiles occur primarily in the outer parts of cluster profiles, where our signal-to-noise is 
low. The similarity between model profiles also means that, in general, the best-fit models of the 
two families have very similar x 2 , with no strong systematic preference for one model or the other. 
Typical x 2 values are 85-90; with ~ 20 datapoints and 3 or 4 model degrees of freedom, the resulting 
reduced values are x 2 ~ 6. This indicates that the uncertainties produced by integrated photometry 
are likely underestimates, and one reason may be that these uncertainties do not account for the 
uncertainty in the background level. Rather than modify the uncert ainties to achieve y 2 ~ 1, we 
modified our use of x 2 i n computing parameter uncertainties (see also lMcLaughlin et al.ll2008l ). We 
scaled the reduced x 2 values such that the best-fit model had x 2 = 1- The 68% confidence limits 
on the parameters are then the minimum and maximum values found in the set of models with 
X 2 < 2. This rescaling gives more realistic estimates of the parameter uncertainties than would 
otherwise be the case. 



How robust are the physical parameters derived from our model fits? One way to estimate this 
is to compare various fits to the same cluster. Although Wo and ro have slightly different meanings 
in King and Wilson models and cannot be directly compared, some derived quantities such as the 
half-light radius and total luminosity are directly comparable. For all clusters we have profile data 
in two different bandpasses, although the outer parts of the profile, derived from number counts, 
are the same in both. There are physical reasons why profiles might change with wavelength (e.g., 
mass segregation, differential reddening), but comparison of model fits in different filters is a useful 
sanity check. Figure [2] shows this comparison: the sca tter between fi l ters is 0.2-0.3 dex. A similar 
comparison between fits for M31 globular clusters by iBarmbv et al.l (|2007l ) found a much smaller 
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scatter, probably because that work analyzed bright clusters, using much deeper data. Figure [2] 
also compares Rh and Ly between Wilson and King models. The scatter is again rather large, 0.15- 
0.25 dex, with the Wilson models offset to larger values. To some extent this is to be expected, 
since Wilson models have larger halos; however some of the Wilson model values (e.g., Rh > 50 pc 
for B015D, B257D, B321, B376, and B448) are physically implausible, because the model-fitting 
resulted in a very large values of the central potential Wo- We do not completely understand the 
reason for this but speculate that it may be related to the combination of the additional power in 
the haloes of Wilson models and the low signal-to-noise of the profiles in the same region. These 
results indicate the limitations of our relatively shallow data, and the limited precision of the model 
measurements will need to be kept in mind during the following analysis. 

For the analysis in the remainder of this paper, we use only a single set of model parameters 
per cluster. Because the King models have fewer implausible values, and also somewhat less scatter 
between filters, we use on the King model parameters for the present cluster sample. Our results 
in Appendix [A] indicate that King model fits may be more robust than Wilson model fits in the 
case where background levels are uncertain, even where the underlying cluster profile is actually a 
Wilson model. Using King models al so allows us to compare th e present sample to the combined 



sample of M31 globulars analysed in iBarmbv et al.l (|2002l . 120071 ): all of that sample has King fits 



while only about one third has Wilson-model fits. Because the focus of this paper is the young M31 
clusters, dominated by blue stars, we use the F439W or F450W-band measurements in preference 
to those from the redder filters. 

The left panel of Figure [3] shows the properties of the present sample of clusters as a function of 
luminosity. Four clusters (vdBO, B327, B342, B368) stand out as having very high central surface 
brightnesses; all except B327 also have correspondingly high concentrations. Figure Q] shows that 
the cores of these clusters do not appear to be resolved in our data. This could be due to the short 
exposure times: if the central cluster light is dominated by a few bright stars, the true integrated 
profile could be very difficult to recover. Structural parameters for these clusters are uncertain. 
Figure [T] also shows that the three M31 young clusters with the largest inferred half-light radii 
(B015D, B321, B448) have relatively low contrast against the resolved stellar background of M31, 
so it is possible that the number counts include some field stars and the resulting Rh values are 
overestimates. The old cluster NB16 has a much smaller Rh and total luminosity than the other 
members of the sample: this cluster is projected on the M31 bulge and its outer stars may be lost 
against the bright background. These issues highlight the limitations of our dataset for the kind 
of structural analysis we are attempting, but the generally good match of model profiles with the 
observational ones gives us confidence that the cluster parameters we measure are reasonable. 

Analyzing the physical properties of M31 young clusters requires converting the observed 
flux-based measurements to luminosities and mass- linked quantities. Conversion from luminos- 
ity to mass is don e using l/-band ma ss-to- light ratios from the population synthesis models of 



Bruzual Chariot] (120031 ) , assuming a IChabrierl (J2003J) IMF and solar meta llicity for all but the 



oldest clusters. Table [T] lists the assumed ages for all clusters: those given bv lPerina et al.l (|2009bl ) 



- 9 - 



for the young clusters from GO-10818, by I Williams &c Hodgd (|200ll ) for the clusters from GO-8296, 
and assumed ages of 13 Gyr for the clusters B083, B222 and B347, B374 and NB16. We assume 
uncertainties of 10% in M/Ly and propagate these through the parameter estimates. While using 
M/Ly ratios determined directly from measured velocity dispersions would avoid the reliance on 
models, velocity dispersions are not available for most of the M31 clusters considered here. The 
use of a single set of population synthesis models also facilitates comparison of clusters in different 



galaxies; the comparison data for other galaxies, (IMcLaughlin fc van der Marelll2005l : lBarmby et al 



2007; 



McLaughlin et alJl2008l ) also used the same model mass-to-light ratios. Tables S] and [5] give 



various derived par ameters for the best-fitt ing models for each cluster (the detail s of their calcu- 
lation are given by IMcLaughlin et al.l 120081 ) . Recently, IKruiissen fc Lamersl (j2008l ) have discussed 
of star cluster mass-to-light ratios due to preferential loss of low-mass stars w ith cluster age. This 
effect is expected to be most important for old clusters, and we have used the IKruiissen &; Lamers 
models to confirm that the change in M/L for young clusters is minimal (< 20%). Since our focus 
in this paper is the young M31 clusters, we therefore do not correct for this effect. 



3. Discussion: young and old clusters in M31 and other galaxies 



Using star clusters as markers of the history of galaxies is aided by knowing how the clusters' 
structural properties change with age and environment. Although absolute ages of star clusters 
are notoriously difficult to determine, relative ages a re more straightforward, and all of the cluster s 
in our sample have ages estimated by CMD fitting ([Williams &; Hodgd l200ll ; iPerina et al.ll2009bl ). 
Can we see evidence for changes in cluster properties with age? In the right panel of Figure [3j 
structural properties for the M31 young clusters are shown as a function of estimated age. None of 
the properties plotted depends on mass-to-light ratio, which is strongly dependent on age. Although 
our sample is small and covers a limited range in age, there is an interesting hint that central surface 
brightness becomes fainter and concentration decreases as a ge increases. This i s cons istent with 
the increase in core radius with age for MC clusters noted by lMackey fe Gilmord (120031 ). Figured] 
explores this further by plotting hq, c, R c , and central mass density po for both the M31 young 
clusters and young clusters in the Magellanic Clouds. While the MC clusters also show a trend for 
central surface brightness to fade with age, it is much weaker than the trend implied by the M31 
clusters alone, and the high-surface-brightness M31 clusters appear to be outliers (possibly artifacts 
due to the limited spatial resolution). Since the central mass density shows very little trend with 
age, the central surface brightness trend is likely due to fading of stellar population and the (weak) 
increase of core radius with age. The dashed line in t he central surface brightne ss panel shows 
the effects of m ass-to-light ratio change predicted by the lBruzual fc Charlotl (J2003J) models with a 
Chabrierl (J2003J) IMF and solar metallicity; the slope shows a reasonable match to the cluster trend. 

Figure d] shows that, with a few exceptions, the young M31 clusters have similar spatial struc- 
ture to young clusters in the Mage llanic Clouds. A number of young massive clusters have recently 
been identified in the Milky Way; iPfalznerl (|2009i ) compiled size and mass measurements of these 
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clusters (|Figerll2008l ; IWolff et al.ll2007l ) to argue that cluster evolution occurs along two well-defined 
tracks in the density-radius plane. Using the conversion between Milky Way cluster size measure- 
ments and half-light radii described in Appendix [A] we have compared cluster half-light radii and 
ages for the young Milky Way clusters together with the M31 and MC clusters in Figure [5j The 
M31 and MC clusters have similar sizes to the 'leaky' Milky Way clusters but lie on the extrapola- 
tion of the age-Rh trend of the 'starburst' MW clusters. This suggests that the starburst clusters 
(which tend to be more massive) are perhaps closer to being analo gs of the young massive clusters 
in other galaxies. We speculate that the two evolutionary paths of iPfalznerl (|2009l ) may be simply 
due to extinction effects, with the 'starburst' clusters having left their host cocoon and the 'leaky' 
clusters still affected by excessive extinction in their outer regions (projection effects may also be 
important). This would imply that starburst clusters are more easily identified in external galax- 
ies, explaining the reasonable match between extragalactic young clusters and Milky Way starburst 
clusters. 

An important question in the study of young massive clusters is whether they will eventually 
become old massive clusters resembling the globular clusters we see today in the Galaxy. Once 
formed, star clusters have no easy way to gain mass , but they do h a ve a number of ways t o lose 



mass or even be completely disrupted (jSpitzer 



1987 



Vesper ini 



1998 



Lamers & Gieled I2006T ) . We 



have computed dissolution times for our cluster sample considering the effects of both the stellar 
and dynamical evolution of star clusters through time. These calculations explicitly account for age, 
metallicity, and half-light radius of all sample star clusters, and treat the effects of evaporation of 
low-mass stars, mass loss due to stellar evol ution, encounters wi th spi ral arms and gian t mole cular 



clouds following in part the prescriptions of lLamers et al.l (|2005l ) and lLamers &: Gieled (120061 ) The 



results are shown in Figured All clusters have dissolution time greater than their ages; however, for 
2 young clusters (B321, B342) and the old cluster B374 these quantities are nearly equal, suggesting 
that they are in the process of dissolving. On average, the young clusters' dissolution times are 
too short to expect them to become old (> 10 10 yr) clusters. However, a few have td > 1 Gyr and, 
if they avoid collisions with giant molecular clouds, might survive to become sparse old globulars. 



Gnedin & Ostriker 


1997; 


Gieles 


2009 


). 


those discovered by 


Krienke & Hodee 



and lHodge et al.l (J2009J), would be even more likely to dissolve. 



Work to date suggests that the structural parame ters of old star clus ters in several nearby 
galaxies show only a weak dependence on environment (jBarmbv et al.l 120071 ). and the comparisons 
above indicate that young clusters in different galaxies are also similar. How do young and old 
clusters compare? Figure [7] shows cluster properties as a fun ction of mass for M31 young cluster s, 
Magellanic Clou d young clusters and Milk y Way globulars ([McLaughlin van der Marell 120051 ) , 



M31 globulars (jBarmbv et al 



2002 



20071 ). and recently-discovered extended M31 halo clusters 
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(|Huxor et alJbood lFI The joint mass- age distribution of the clusters differs by galaxy: some of this 
is due to complex selection effects (e.g., the M31 globular sample is incomplete and biased toward 
more luminous clusters, and the sample of Milky Way YMCs i s also incomplete), but there are hints 
of real differences between galaxies; see iPerina et al.l (|2009bl ) for a more detailed discussion. The 
properties of the five old clusters in our sample are similar to those of M31 and Milky Way globulars, 
while the properties of M31 young clusters overlap with those of both the young Magellanic Cloud 
clusters and the low-mass Milky Way globular clusters. Thus the M31 young clusters do not appear 
to be fundamentally different types of object from those already known. On average, the younger 
clusters have larger sizes and higher concentrations (where larger c implies a larger tidal radius 
for the same scale radius) than old clusters of the same mass. The young clusters therefore have 
larger tidal radii, which makes them more susceptible to dynamical destruction: small-r^ clusters 
are more likely to survive to old age. The larger spread in properties of low-mass clusters compared 
to higher- mass clusters may indicate lower data quality for these fainter objects, rather than an 
intrinsic difference in properties. 

By now it is well-known that old star clusters i n the Milky Way and other galaxies des cribe 
a 'fundamental plane' (FP) in structural properties (|Djorgovskilll995l ; iDiorgovski et al.lll997l ). al- 
though the separatio n of clusters from oth er t ypes of objects has become less well-defined in recent 
years. The results of lBastian et al.l (|2006l ) and iKissler-Patig et al.l (|2006l ) indicate that young mas- 
sive clusters fall on similar fundamental planes to those of old clusters. Those results make use of 
cluster velocity dispersions, while in this work, we must use mass-to-light ratios from population 
synthesis models applied to the photometry instead of inde pendent mass estima tes. The upper-right 
panel of Figure [7] shows one view of the FP, as defined by iMcLaughlinl ([2000). The old clusters in 
our sample fall nicely on this relation, as do most of the younger clusters. The observed correlation 
between mass and binding energy Ef, is expected, since by definition = f{c)M 2 /Rh where /(c) 
is a weak function of cluster concentration c. However, the tightness of the correlation shows that 
there is very little relation between young cluster mass and (see also lower-right panel), and no 
offsets in the basic properties of the cluster shapes between old and young clusters. 

Figure shows a different view o f the fundamental plane, more akin to t he parameters usually 
shown for elliptical galaxies (see also iMcLaughlin 2003; IStrader et al.ll2009l ). The left two panels 
show the surface-brightness-based fundamental plane relations, with a large offset between the 
young M31 and MC clusters (light grey symbols) and the old clusters. This is to be expected 
because of the young clusters' lower mass-to-light ratios. When we instead plot quantities related to 
the mass density (right panels), the young clusters fall on the same relations as the old clusters. The 
tightness of the relations primarily reflects the use of mass-to-light ratios to compute both central 
velocity dispersion ctq and mass density S. Again, however, the lack of offset and similar scatter 



4 Mass measurements for all clusters are derived using mass-to-light ratios. As discussed in £12.41 these ratios are 
affected by cluster dyn amical evo l ution. Correcting for this effect is non-trivial and beyond the scope of this paper; 
however the results of Kruiissen ( 20081 ) imply that doing so would increase the spread of the old clusters' mass 
distribution and shift it to lower masses. 
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between the young and old clusters confirms t heir similar overall s tructures. Recent measurements 
of M31 globular clusters' mass-to-light ratios (jStrader et al.ll2009l ) have shown that these clusters 
do follow the FP relations as expected from model mass-to-light ratios. Similar measurements for 
young clusters should show whether young clusters do the same. If so, this would indicate that the 
FP reflects conditions of cluster formation and is not merely the end product of cluster dynamical 
evolution. 



Bonatto Bical (|2005l ) argue that Milky Way open clusters fall on a plane in the three- 



dimensional space of total mass, core radius, and projected core mass density. We can compare 
this space to the FP using with an approximate relation between mass and central velocity disper- 
sion. The least-squares fit for the young MC clusters (the most populous sampl e of yo ung clusters 
available) gives logoo = 0.34 log M — 1.38; combined with the iBonatto k, Bical (120051 ) cluster pa- 
rameters, we find that the Milky Way open clusters fall approximately on the other young clusters 



in the top right panel of F igure El This suggests that the Milky Way open 



with S ~ 10 2 M pc~ 

cluster plane indicated by IBonatto &: Bical (120051 ) may in fact be the same FP defined by other 
star clusters, which ha ve projected mass densities higher by up to four orders of magnitude. As 
Bonatto &; Bical (|2005l ) discuss, this result remains to be confirmed with large samples, but it is 
certainly intriguing in its implications for a 'universal' star cluster fundamental plane. 



4. Summary and Directions for Future Work 

This series of papers has established that a sample of candidate young star clusters in M31 are 
indeed young, massive clusters, with properties similar to those of other young clusters in Local 
Group galaxies. Our current data does not allow us to detect the extended haloes characterized by 
Wilson models and seen in other young clusters; the more compact King models provide adequate 
fits to the data. The structural parameters measured in this paper show the M31 clusters to be 
typical young clusters, with masses of 10 4-5 M Q , half-light radii of 3-20 pc, and dissolution times 
of < 5 Gyr. While the basic similarity between young clusters in different Local Group galaxies, 
and between young and old clusters, seem well-established, many questions remain. What is the 
precise form of the age-size relation? Do cluster mass-to-light ratios evolve with age as predicted by 
dynamical and stellar evolution models? What fraction of the stellar disk in galaxies is comprised 
of dissolving clusters? Is there a relation between the cluster formation and local star formation 
rate, or other galaxy properties? Large cluster samples with high-quality data will be needed to 
address these and other questions about the relationship and history of star clusters and their 
parent galaxies. 

Facilities: HST (WFPC2) 



We thank D.E. McLaughlin for the use of the GRIDFIT software. P.B. acknowledges research 
support through a Discovery Grant from the Natural Sciences and Engineering Research Council of 



-13- 



Canada and an Ontario Early Researcher Award. S.P. and M.B. acknowledge the financial support 
of INAF through PRIN 2007 grant no. CRA 1.06.10.04 "The local route to galaxy formation." 
J.G.C. acknowledges support from NASA from grant GO-10818. T.H.P acknowledges support 
through the Plaskett Research Fellowship at the Herzberg Institute of Astrophysics. 



Artificial cluster tests 



Deriving surface brightness profiles of star clusters in Local Group galaxies requires careful 
analysis. The clusters are only partially resolved into individual stars, and they are observed 
together with a galactic background which may also be resolved into stars. The purpose of this 
section is to investigate the best methods for extracting structural parameters of 'semi-resolved' 
clusters, particularly from relatively shallow images, and to quantify the uncertainties of those 
parameters. This can best be done by analyzing profiles de rived from images of artific ial clusters 
whose structural parameters are known. A related study by iNovola Gebhardtl (|2006l ) simulated 
integrated photometry from HST observations of Galactic GCs; however th e focus of that study wa s 
on recovering the structure of cluster cores rather th an overall s tructure. iBonatto &: Bical ([2008) 
also carried out a similar study, but considering only iKingl (|1962l ) models for Galactic clusters. 



The first step in analyzing simulated star cluster profiles is to determine the type of model 
profile and range of parameter space to be covered. The analysis of iMcLaughlin &: van der Marel 
(|2005l ) showed that Wilson models were adequate to describe both Milky Way and Magellanic 
Cloud cluster profiles, so we chose this set of models for our artificial clusters. Since we are 
interested in differences between young and old clusters we examined the distribution of scale 
radius rn and central potential Wo for both young and old Magellanic Cloud clusters as given by 
McLaughlin &: van der Marell (120051 ): Wq ranged from 1 to 10 with a typical value W$ ~ 5 while 
ro ranged from 0.2 to 20 pc with a typical value ro ~ 2 pc. The range of implied half-light radii is 
1-35 pc. 

Our artificial clusters were generated from Wilson profiles with 8 values of ro between 0.5 and 
11 pc, and 9 values of Wq between 2 and 10. For each {Wq,tq) pair we generated clusters with four 
different population sizes: iV* = 100, 300, 1000, 3000. The stars' projected spatial positions were 
generated by selecting the projected radial coordinate from the probability distribution associated 
with the Wilson profile 



p{R) 



(Al) 



^(R')R'dR' 

and generating the angular coordinate 9 at random. The stars' luminosities were generated by 
selecting from an observed 'young cluster' luminosity distribution, uncorrected for completeness. 
The distribution was generated by combining the observed magnitudes of stars in the four most 
populous clusters in the GO-10818 program (VdBO, B257D, B475, B327). Separate luminosity 
distributions were used in each of the two observational bands. 
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The specific observations being modeled are the same as those in the GO-10818 program. We 
generated images of the simulated clusters by inserting artificial stars modeled with the appropriate 
PSF near the center of a WFPC2/PC image of a field in M31. The background images used were 
the observations of 'B195D' from the GO-10818 program; the PC chip w as essentially empty in this 



observation because of an error in the input coordinates (for details, see lPerina et al.ll2009al ). This 
field is located in the south-west disk of M31. Figure shows a sample of the simulated cluster 
images, together with some sample M31 clusters for comparison. The simulated clusters cover a 
wider range of properties than the real clusters: some of the simulated clusters were in fact not 
visually apparent in the images. These 'clusters' had few stars (iV* = 100 or N* = 300) and very 
large half-light radii, more akin to dwarf galaxies than to objects recognizable as star clusters. They 
are not considered further in this analysis. 

Surface density profiles for the simulated clusters were derived in several different ways. The 
first method ('number counts'), derived the surface density as simply the number of stars per unit 
area in annular bins. Since the locations of all stars are known precisely for the simulated clusters, 
this method represents the best possible data for surface density profiles. Deriving structural 
parameters from such data tests the fitting routine itself and also the extent to which density 
profiles can be derived from a limited number of stars. Stars were counted in overlapping annular 
bins of width 3 pixels (0.5 pc) inside a radius of 20 pixels (3.4 pc) and width 10 pixels (1.7 pc) 
outside this radius. 

For real star clusters, crowding limits the ability to resolve individual stars and hence derive 
surface density profiles through number counts. We also derived surface density profiles of clusters 
using isophota l phot ometry with the IRAF ellipse package, similar to the method described in 



Barmbv et al.l (|2007l ). We refer to this as the 'integrated photometry' method. We also combined 



the number count an d integrated photometry methods in a 'hybrid' method similar to that used by 



Federici et al.l (|2007l ) . This involves matching the intensity scales of the two profiles by fitting both 
profiles to smooth curves in the region r = 5 — 10 pc. The switch-over from integrated photometry 
to number counts was made at a radius of 7 pc (40.6 pixels), where in general both types of profile 
had good signal to noise. 

Wilson models were fit to the artificial cluster data using the GRIDFIT program described 
in §2.31 As for the real clusters, instrumental PSF profiles were convolved with the model profiles 
before comparison to the data. Unlike the real clusters, however, the background level for the 
artificial clusters was fixed at zero. For clusters of all sizes, the number count input returned 
fitted parameters in good agreement with the input parameters. The offsets between input and 
output parameters are (mean ± standard error) A Wo = (Wo,i n — Wo )0 ut)/Wo,in = 0-06 ± 0.02 
and Aro = (ro,m — ^o,out)/^*o,in = —0.13 ± 0.03 pc. As expected, the larger- iV* clusters return 
more accurate values, with scatter 2-3 times lower for iV* = 3000 than for N* = 300 clusters. 
Figure [10] compares the best-fit and input structural parameters of the simulated clusters for the 
integrated photometry and hybrid methods. Particularly for clusters with larger input ro, integrated 
photometry alone tends to result in overly-large values of Wq and overly-small values of ro- For 
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these clusters, the distinction between profiles of different Wq occurs at a point in the radial profile 
where the density of stars is too low for the ELLIPSE algorithm to converge. The addition of number 
count data beyond this point improves the fit, as the figure shows. For integrated photometry alone, 
AWo = -0.56 ± 0.07 and Ar = 0.24 ± 0.04 pc; for the hybrid method, AW = -0.02 ± 0.02 and 
Ar = -0.05 ± 0.03 pc. 

When fitting model profiles to cluster data, the correct model family is not not known a 
priori. What happens if artificial 'Wilson' clusters are fit with King models instead? We tried 
this experiment with our artificial clusters and were surprised to find that, except for a handful of 
objects, the two model families returned nearly identical \ 2 values: the median fractional difference 
(Xk ~ Xw)/^w = 0-01- While the meaning of model parameters such as the scale radius ro differs 
between model families, some derived quantities such as the core and half-light radii (R c ,Rh- see 
Table |4] for description) are directly comparable. Figure [IT] shows this comparison. There is 
very good agreement between the two model families in measurements of core radii, and reasonable 
agreement in measurements of half-light radii. The agreement in Rh is poorer for the largest clusters 
(Rh > 20 pc, a larger size than usually seen in real clusters), where t he King models return smaller 



sizes than the Wilson models. This is consistent with the results of iMcLaughlin &: van der Mare! 



hood ) who found that the two model families gave generally consistent results for Milky Way and 



Magellanic Cloud clusters as long as the radius of the last data point -Ri as t ^ 5Rh- 

The situation of observational profiles with a limited radial range bears further investigation. 
The analysis of simulated clusters to this point has not considered the effects of background level 
fluctuations. The GRIDFIT code is able to fit a constant background level added to the intensity 
profile, and we verified through simple experiments that input values were correctly recovered. 
However, the limitations of short exposures and small-number statistics suggest that determining 
the correct background level — and thus being able to correctly trace cluster profiles out to large 
projected radii — will be much more difficult for the real cluster data. We therefore experimented 
with removing points in the profile data beyond -Ri as t = 1> 2, and 5Rh (where Rh was computed 
from the input model profile) and fitting both King and Wilson models to the remaining points. As 
expected, recovery of the input cluster parameters was better for the more extensive profiles, for 
both model families. For i?i as t — 1 ; both model families returned Rh values that were, on average, 
larger than the input. Some model fits were 'catastrophic failures', with Rh{out) > 2i?/j(in); 
this situation usually occurred for clusters where the number of profile data points was < 10. 
Interestingly, for all three values of i?i as t, King model fits had fewer catastrophic failures than 
Wilson models, and also slightly smaller scatter in the difference between fit and true parameters. 
Since the primary difference between King and Wilson model profiles is the more extended halo of 
the latter, this suggests that King models may be a better choice for fitting noisy cluster profiles. 

Finally, we considered the issue of comparison between different measurements of star cluster 
size. While Milky Way globulars and extragalactic clusters are most often characteriz ed with half 



light or core radi i, recent complilations of data for massive young Milky Way clusters (jFigerl 120081 : 



Wolff et al.ll2007l ) measure cluster size as the mean or median distance ((R) or R) of the cluster stars 
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from the geometric centroid. Since these young Milky Way clusters may well not be dynamically 
relaxed ([Goodwin &: Bastianl 120061 ) , it may not make sense to fit the same types of dynamical 
models to them as to old clusters, but it is still desirable to find a way to compare sizes between 
groups of clusters. Since we know the positions of all stars in our artificial clusters, we can easily 
compute the statistical measurements of size for our model clusters, and compare them to (model 
values of) R c and Rh- (R) and R are very well-correlated for all of our model clusters, with a 
best-fit linear relation R = 0.67(i?) — 0.36. The correlation between (R) and R c is rather poor 
(unsurprising as R c depends critically on the exact shape of the cluster profile), but there is a good 
match between (R) and Rh for models which are not too extended (Wq < 6). Figure [T21 shows the 
data and least-squares fits: (R) = 0.77Rh + 0.23, and R = 0.53Rh + 0.10. We conclude that, with 
some scaling, the mean or median projected separation of stars from a cluster center are reasonable 
proxies for the half-light radius. 



REFERENCES 

Barmby, P., Huchra, J. P., Brodie, J. P., Forbes, D. A., Schroder, L. L., & Grillmair, C. J. 2000, 
AJ, 119, 727 

Barmby, P., Holland, S., & Huchra, J. P. 2002, AJ, 123, 1937 

Barmby, P., McLaughlin, D. E., Harris, W. E., Harris, G. L. H., & Forbes, D. A. 2007, AJ, 133, 
2764 

Bastian, N., Saglia, R. P., Goudfrooij, P., Kissler-Patig, M., Maraston, C., Schweizer, F., & Zoccali, 
M. 2006, A&A, 448, 881 

Bendinelli, O., Cacciari, C., Djorgovski, S., Federici, L., Ferraro, F. R., Fusi Pecci, F., Parmeggiani, 
G., Weir, N., & Zavatti, F. 1993, ApJ, 409, L17 

Bica, E., Bonatto, C., Dutra, C. M., & Santos, J. F. C. 2008, MNRAS, 389, 678 

Bonatto, C. & Bica, E. 2005, A&A, 437, 483 

Bonatto, C. & Bica, E. 2008, A&A, 477, 829 

Brodie, J. P., Schroder, L. L., Huchra, J. P., Phillips, A. C., Kissler-Patig, M., & Forbes, D. A. 
1998, AJ, 116, 691 

Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 

Caldwell, N., Harding, P., Morrison, H., Rose, J. A., Schiavon, R., & Kriessler, J. 2009, AJ, 137, 
94 

Chabrier, G. 2003, PASP, 115, 763 



-17- 



Cohen, J. G. & Freeman, K. C. 1991, AJ, 101, 483 

Cohen, J. G., Matthews, K., & Cameron, P. B. 2005, ApJ, 634, L45 

Djorgovski, S. 1995, ApJ, 438, L29 

Djorgovski, S. & King, I. R. 1986, ApJ, 305, L61 

Djorgovski, S. & Meylan, G. 1994, AJ, 108, 1292 

Djorgovski, S. G., Gal, R. R., McCarthy, J. K., Cohen, J. G., de Carvalho, R. R., Meylan, G., 
Bendinelli, O., & Parmeggiani, G. 1997, ApJ, 474, L19 

Elson, R. A. W., Fall, S. M., & Freeman, K. C. 1987, ApJ, 323, 54 

Federici, L., Bellazzini, M., Galleti, S., Fusi Pecci, F., Buzzoni, A., & Parmeggiani, G. 2007, A&A, 
473, 429 

Figer, D. F. 2008, in IAU Symposium, Vol. 250, IAU Symposium, ed. F. Bresolin, P. A. Crowther, 
& J. Puis, 247 

Fusi Pecci, F., et al. 1994, A&A, 284, 349 

Fusi Pecci, F., Bellazzini, M., Buzzoni, A., De Simone, E., Federici, L., & Galleti, S. 2005, AJ, 130, 
554 

Galleti, S., Bellazzini, M., Federici, L., Buzzoni, A., & Fusi Pecci, F. 2007, A&A, 471, 127 

Galleti, S., Federici, L., Bellazzini, M., Fusi Pecci, F., & Macrina, S. 2004, A&A, 416, 917 

Gieles, M. 2009, MNRAS, 394, 2113 

Gnedin, O. Y. & Ostriker, J. P. 1997, ApJ, 474, 223 

Goodwin, S. P. & Bastian, N. 2006, MNRAS, 373, 752 

Grillmair, C. J., Ajhar, E. A., Faber, S. M., Baum, W. A., Holtzman, J. A., Lauer, T. R., Lynds, 
C. R., & O'Neil, Jr., E. J. 1996, AJ, 111, 2293 

Hasegan, M. et al. 2005, ApJ, 627, 203 

Hodge, P. W., Krienke, O. K., Bellazzini, M., Perina, S., Barmby, P., Cohen, J. C, Puzia, T. H., 
& Strader, J. 2009, AJ, in press (|arXiv:0907.1305p 

Holland, S., Fahlman, G. G., & Richer, H. B. 1997, AJ, 114, 1488 

Holtzman, J. A., et al. 1992, AJ, 103, 691 

Hubble, E. 1932, ApJ, 76, 44 



-18- 



Huxor, A. P., Tanvir, N. R., Irwin, M. J., Ibata, R., Collett, J. L., Ferguson, A. M. N., Bridges, 
T., & Lewis, G. F. 2005, MNRAS, 360, 1007 

Jordan, A. 2004, ApJ, 613, L117 

Kharchenko, N. V., Piskunov, A. E., R6ser, S., Schilbach, E., & Scholz, R.-D. 2005, A&A, 438, 
1163 

King, I. 1962, AJ, 67, 471 
King, I. R. 1966, AJ, 71, 64 

Kissler-Patig, M., Jordan, A., & Bastian, N. 2006, A&A, 448, 1031 

Krienke, O. K. & Hodge, P. W. 2007, PASP, 119, 7 

-. 2008, PASP, 120, 1 

Kruijssen, J. M. D. 2008, A&A, 486, L21 

Kruijssen, J. M. D. & Lamers, H. J. G. L. M. 2008, A&A, 490, 151 
Kundu, A. 2008, AJ, 136, 1013 

Lamers, H. J. G. L. M.„ Gieles, M. , Bastian, N., Baumgardt, H., Kharchenko, N. V. & Portegies 
Zwart, S. 2005, A&A, 441, 117 

Lamers, H. J. G. L. M. & Gieles, M. 2006, A&A, 455, L17 

Larsen, S. S. 2004, A&A, 416, 537 

Larsen, S. S. & Richtler, T. 1999, A&A, 345, 59 

McLaughlin, D. E. 2000, ApJ, 539, 618 

Mackey, A. D., & Gilmore, G. F. 2003, MNRAS, 338, 85 

McLaughlin, D. E. 2003, in Astronomical Society of the Pacific Conference Series, Vol. 296, New 
Horizons in Globular Cluster Astronomy, ed. G. Piotto, G. Meylan, S. G. Djorgovski, & 
M. Riello, 101 

McLaughlin, D. E., Barmby, P., Harris, W. E., Forbes, D. A., & Harris, G. L. H. 2008, MNRAS, 
384, 563 

McLaughlin, D. E. & van der Marel, R. P. 2005, ApJS, 161, 304 
Meylan, G. & Djorgovski, S. 1987, ApJ, 322, L91 
Noyola, E. & Gebhardt, K. 2006, AJ, 132, 447 



-19- 



Perina, S., et al. 2009a, A&A, 494, 933 
Perina, S. et al. 2009b, in prep 
Pfalzner, S. 2009, A&A, 498, L37 

Piskunov, A. E., Schilbach, E., Kharchenko, N. V., R6ser, S., & Scholz, R.-D. 2007, A&A, 468, 151 

-. 2008, A&A, 477, 165 

Racine, R. 1991, AJ, 101, 865 

Racine, R. & Harris, W. E. 1992, AJ, 104, 1068 

Rich, R. M., Mighell, K. J., Freedman, W. L., & Neill, J. D. 1996, AJ, 111, 768 

Schilbach, E., Kharchenko, N. V., Piskunov, A. E., Roser, S., & Scholz, R.-D. 2006, A&A, 456, 523 

Sersic, J.-L. 1968, Atlas de Galaxias Australes (Cordoba: Observatoria Astronomico) 

Spitzer, L. 1987, Dynamical evolution of globular clusters (Princeton, NJ, Princeton University 
Press) 

Stanek, K. Z. & Garnavich, P. M. 1998, ApJ, 503, L131 

Strader, J., Smith, G. H., Larsen, S., Brodie, J. P., & Huchra, J. P. 2009, AJ, 138, 547 
van den Bergh, S. 1994, AJ, 108, 2145 
Vesperini, E. 1998, MNRAS, 299, 1019 

Weidner, C, Kroupa, P., & Larsen, S. S. 2004, MNRAS, 350, 1503 
Whitmore, B. C. & Schweizer, F. 1995, AJ, 109, 960 
Williams, B. F. & Hodge, P. W. 2001, ApJ, 548, 190 
Wilson, C. P. 1975, AJ, 80, 175 

Wolff, S. C, Strom, S. E., Dror, D., & Venn, K. 2007, AJ, 133, 1092 



This preprint was prepared with the AAS IATgjX macros v5.2. 



-20- 



0.1 1 10 100 



R/pc 

0.1 1 10 100 



10 4 
10 3 
10 2 
10 

10 4 
10 3 
10 2 

i io 



< 10 4 
10 3 
10 2 
10 

10 4 
10 3 
10 2 
10 




- B043 n 

I I I L 



B081 



I I I L 



I I I L 



I I I L 



+- B066 

I I I L 



singes 



ii- B083* 



v 



I I I L 




r B222* ir -±- B257D -±- 

t I I I i~tt I I I inF I I I i -Fr I I I ud 



I I I 



ffi m fTn ID ffi 
— m x >— i 



I I I L 



0.1 1 10 0.1 1 10 0.1 1 10 0.1 1 10 

R/arcsec 



Fig. 1. — M31 cluster surface brightness profiles together with the best-fitting models. Each cluster 
is shown in two sub-panels, with the bluer filter (F439W or F450W) on the left and the redder 
filter (F555W or F814W) on the right . Clus ters with an asterisk after their names are likely to be 
old. Black lines are best-fitting iKina (|1966l ) models; grey lines (most are directly over the black 
lines) are best-fitting IWilsonl ()1975l ) models. Solid lines are model profiles after convolution with 
the PSF; dash-dot lines are profiles before convolution. 
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Fig. 1. — Continued 
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Continued. Note that the last four clusters are plotted with a different vertical scale. 
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Fig. 2. — Comparison of half-light radii and total luminosity (converted to the V band) for Wilson 
and King models fit to surface brightness profiles of M31 young clusers. Bottom: comparison 
between observations of the same cluster in different filters (hexagons: King models, stars: Wilson 
models). Top: Comparison of Wilson and King model fits to the same cluster (squares: red filter, 
triangles: blue filter). 
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Fig. 3. — Concentration index, central surface brightness, and half-light radius for M31 young 
clusters as functions of total model luminosity (left) and estimated age (right). The old clusters 
are shown with gray symbols; although their ages are estimated at > 10 10 yr, they are plotted at 
10 9 yr in the right panel to condense the horizontal axis scale. 
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Fig. 4. — Concentration index, central surface brightness, and central mass density for M31 
(squares) and Magellanic Cloud (triangles) young clusters as functions of estimated age. The 
dashed line in the central surface brightness panel shows the expected change in surface brightness 
due to changes in mass-to-light ratio with age (vertical normalization is arbitrary). 
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Fig. 5. — Young star cluster ages and sizes. Squares: M31 clusters from the present sa mple 
triang les: young Magel l anic C loud clusters; circles: young massive Milky Way clu sters from iFiger 
(l2008h and lWolff etail ^Ooj). The two groups of Milky Way clusters identified bv lPfalzneri fcood ) 
are labeled. 
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Fig. 6. — Dissolution times for M31 star clusters, compared to cluster ages. Four of the five old 
clusters are plotted at the same position, with dissolution times 20 Gyr and assumed ages 13 Gyr. 
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Fig. 7. — Structural properties of young and old star clusters in M31, young clusters in the Mag- 
ellanic Clouds, and globular clusters in the Milky Way, shown as a function of cluster mass. Top 
left: central surface brightness; top right: binding energy; lower left: concentration; lower right: 
half-light radius. Filled squares: M31 clusters fro m the present sample (black: young clusters, grey: 
old clusters); open squares: old M31 clust e rs fro m lBarmby et al.l (120071 . 12002 ): stars: 'extended lu- 
minous clusters' in M31 from lHuxor et aL ( 2005 ): filled triangles: young Magellanic Cloud clusters. 
Error bars show median uncertainties for the young M31 clusters. 
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Fig. 8. — Views of the star cluster fundamental plane, with core parameter relations in the 
bottom panels and half-light parameter relations in the top panels. o"o is predicted central velocity 
dispersion and X represents surface mass density either in the cluster core or at the half-light radius. 
Left panels show surface brightness while right panels show mass surface density. Filled squares: 
M31 clusters from the present sample (light grey: young clusters, dark grey: old clusters); open 
triangles: old Magellanic Cloud and Fornax clusters; open circles: Milky Way globulars; crosses: 
NGC 5128 globulars. Other symbols as in Figure [71 Error bars show median uncertainties for the 
young M31 clusters. 
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Fig. 9. — Top row: four M31 star clusters observed as part of program GO-10818 with HST. 
Left to right: B015D, B081, B222, B475. Second row: simulated clusters with central potential 
W = 6 and scale radius r = 2 pc, with (left to right) = 100, 300, 1000, 3000. Third row: 
simulated clusters with = 3000, central potential Wo = 6 and scale radius (left to right) 
ro = 0.5,1.5,5,11 pc. Fourth row: simulated clusters with = 3000, scale radius ro = 2 pc, 
and (left to right) Wo = 2, 4, 8, 10. All images are 800 s exposures in the F450W filter on the 
WFPC2/PC chip; each sub-image is 13.7 x 13.7 arcsec (51.7 x 51.7 pc at the distance of M31). 
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Fig. 10. — Comparison of input and output structural parameters for simulated star clusters. 
The ouput parameters are derived from fitting Wilson models to surface density profiles derived 
from simulated HST/WFPC2 images of the clusters. Left: profiles measured with integrated 
photometry only; right: profiles measured with integrated photometry and number counts; top: 
difference in central potential AWq = (Wo,i n — Wo,out)/Wo,in; bottom: difference in scale radius 
Ar = (r ,i„ - ro,out)/ro,in; 
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Fig. 11. — Comparison of cluster size measurements for fits of model den sity p r ofiles to artificial 
cluster pr ofiles. Top: co re radius R c ; bottom: half-light radius R^, circles: iKina (|1966l ) model fits; 
triangles: IWilsonl (|1975l ) model fits. 
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Fig. 12. — Comparison of model half-light radius Rh to mean and median projected radius for 
artificial clusters. Circles: mean; triangles: median; filled symbols: models with Wq < 6; open 
symbols: models with Wo > 6. Solid line: least-squares fit to filled circles; dotted line: least- 
squares fit to filled triangles. 
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Table 1. Data for M31 young clusters 



Name 3, 


Dataset 1 


Dataset 2 


Filter 1 


Exposure 1 [s] 


Filter 2 


Exposure 2 [s] 


E(B - V) 


log age [yr] 


B015D 


„rt ■ -i a n r-i ol 
Uypil4U[lZJ 


uypi!4U[34J 


r 4D0W 


800 


T7Q 1 A \ XT 

r ol4 W 


800 


0.65 


7.85 


Dn/lfl 

r>U4U 


uypiUOU[lzJ 


uypiUoU[o4J 


r 40UW 


Qnn 
oUU 


T?Q1 AWT 

r ol4W 


Qnn 
oUU 


n oq 


i nn 

i .yu 


B043 


UypiU2z[lZJ 


uypiUZ2[34J 


r 450 W 


800 


r o!4 W 


800 


0.23 


7.90 


r>Ubo 


uypiz4U[lzJ 


uypiz4U[o4J 


T?A ^n\\T 
r 40U W 


oUU 


T?C1 AWT 


Qnn 


n oq 


I.OO 


DflOl 

rSUol 


uypn (U[izj 


uypil / U[o4J 


T?A ^CtWT 


vi in 
oUU 


T?Q1 AWT 

r ol4W 


Qnn 
oUU 


n Qn 


O 1 c 

o.lo 


B083 


.id^ioxnh ol 

uypizoU[izj 


uypizoU[34J 


r 450 W 


800 


T7qi AWT 

r o!4 W 


800 


0.20 


10.11 


dZZZ 


-1 QnTl Ol 

uypiioU[izj 


uypiloU[o4J 


T?A ^ClWT 

r 40U W 


cnn 
oUU 


T?C1 AWT 

r 514 W 


Qnn 
oUU 


n on 
U.zU 


Q nn 

o.yu 


B257D 


u9pil00[12] 


u9pil00[34] 


F450W 


800 


F814W 


800 


0.30 


7.90 


B315 


u5bj010[12] 


u5bj010[78] 


F439W 


1600 


F555W 


1200 


0.31 


8.00 


B318 


u9pi020[12] 


u9pi020[34] 


F450W 


800 


F814W 


800 


0.17 


7.85 


B319 


u5bj020[12] 


u5bj020[78] 


F439W 


1600 


F555W 


1200 


0.23 


8.00 


B321 


u9pil50[12] 


u9pil50[34] 


F450W 


800 


F814W 


800 


0.25 


8.23 


B327 


u9pi030[12] 


u9pi030[34] 


F450W 


800 


F814W 


800 


0.20 


7.70 


B342 


u5bj030[12] 


u5bj030[78] 


F439W 


1600 


F555W 


1200 


0.20 


8.20 


B347 


u9pi230[12] 


u9pi230[34] 


F450W 


800 


F814W 


800 


0.06 


10.11 


B368 


u5bj040[12] 


u5bj040[78] 


F439W 


1600 


F555W 


1200 


0.20 


7.80 


B374 


u9pi070[12] 


u9pi070[34] 


F450W 


800 


F814W 


800 


0.30 


8.80 


B376 


u9pi080[12] 


u9pi080[34] 


F450W 


800 


F814W 


800 


0.30 


8.00 


B448 


u9pi200[12] 


u9pi200[34] 


F450W 


800 


F814W 


800 


0.35 


7.90 


B475 


u9pi090[12] 


u9pi090[34] 


F450W 


800 


F814W 


800 


0.35 


8.30 


NB16 


u9pil20[12] 


u9pi012[34] 


F450W 


800 


F814W 


800 


0.25 


10.11 


V031 


u9pil30[12] 


u9pil30[34] 


F450W 


800 


F814W 


800 


0.35 


8.45 


VDBO 


u9pi010[12] 


u9pi010[34] 


F450W 


800 


F814W 


800 


0.20 


7.60 



a Naming convention of the Revised Bologna Catalog jGalleti et alj|2004h is used. See that work for coordinates. 



Table 2. Calibration data for WFPC2 imaging 



filter 


zcropoint 


M Q 


Conversion factor a 


F439W 


22.987 


5.55 


45.138 


F450W 


23.996 


5.31 


14.274 


F555W 


24.621 


4.83 


5.163 


F814W 


23.641 


4.14 


6.744 



a Multiplicative conversion between surface bright- 
ness in counts s _1 arcsecc" 2 and intensity in 
L pc" 2 . 



Table 3. Basic Parameters of Fits to Profiles of M31 Young Clusters 



Name 


Filter 


N pts 


Model 


X 2 ■ 

Mm ii 


[L 0P c- 2 ] 


W 


c 


Ho 

[mag arcsec -2 ] 


log r 
[arcsec] 


log r 

[PC] 


B015D 
B015D 
B040 


F450W 
F814W 
F450W 


21 
21 
21 


K66 
W 

K66 
W 

K66 
W 


323.12 
386.35 
231.70 
377.92 
44.18 
50.75 


7.5 
7.5 
12.8 
12.8 

33.18 ± 3.56 
21.84 ±5.10 


10.20t^ 

io.80ii;Jg 

14 40+ 1 ' 40 

14.90±iJ° 
q fin +0.40 
9 - ou -0.30 

q sfl +o.50 


n OQ + 0.18 

z " 3y -0.17 

n no+0.13 

°- ao -0.06 
o 9 o+0.31 
°* zo -0.21 
4 1K+0.39 

2 26+ ' 09 

O qr,+0.02 

■ 3 - sz -o.oo 


i6.iit8:S 

12 61 + ' 48 

±Z.O±_ gg 

12.4718; 5 7 J 

15.4418;?? 
15.4818-.™ 


-0.64018:1?! 

-o.65ol8;l?l 

-1.758±8:|?| 
-l.804_ 2g7 
-0 Qfi7+ 048 

u - yD ' -0.067 
-0 Q31+ - 054 

U.MJ1_ Qgg 


-o.06il8;l?I 

-0.07l!8i?l 

-1.17918;^ 

i 99 c+0.215 
l.zzo_ 2g7 
n q S 7+0.048 

-0 S52+ ' 054 



Note. — Table [3] is available in its entirety in the electronic edition of the Journal. A short extract from it is shown here, for guidance 
regarding its form and content. Column descriptions: Xmin : unreduced \ 2 of best-fitting model; ibkg : model-fit background intensity (values 
without uncertainties indicate clusters for which the background was fixed manually); Wo'- model-fit central potential; c = log(rt/ro): model-fit 
concentration (rt is tidal radius, given in Table [4}; fio'- model-fit central surface brightness; log ro: model-fit scale radius. Uncertainties are 
68% confidence intervals, computed as described in the text. 



Table 4. Derived Structural and Photometric Parameters for M31 Young Clusters 



Name 


Filter 


V color 


Model 


log r t j d 


log Rc 


log Rh 


log (Rh/RJ 


log h 


log io 


log Ly 


Vtot 


log I h 






[mag] 




[PC] 


[PC] 


[PC] 




[L 0iV pc" 2 ] 


[£©,VPC" 3 ] 


[Lq,v] 


[mag] 


[Lq, vpc~ 2 ] 


B015D 
B015D 
B040 


F450W 
F814W 
F450W 


-0.114 ±0.1 
0.457 ±0.1 
-0.029 ±0.1 


K66 
W 

K66 
W 

K66 
W 


9 qo+0.06 
z - oo -0.07 
n on+0.07 

■J-^-o.oo 
2 05 +0 - 03 

ZiUO -0.01 
2 qo+0.10 

9 Q7 +0.05 
z - y '-0.05 


U.UDO_q 11Q 

-0-076±g:iS| 

1 17Q+0.196 
L - L '"-0.279 

-1 224+ ' 215 
1.ZZ^_q 2gg 

-0 SQ':i +0 - 047 

-0 S61+ - 052 


1 346+°- 100 

l.Ci4D_ 12Q 

1 74fi+ (M,el 
L - ' 40 -0.051 

i n«fi+ ' 014 

i , 19 + 0.053 

i.oiz_ 025 
n 85Q+ ' 047 

u.ooo_ 

1 9Q9+0.022 
l.z»z_ 032 


1 411+0-210 
i - 4ii -0.226 
1 091 +0-170 
l.OZl_ 16g 

9 9R4+0- 288 
Z - ZD4 -0.194 
9 CO7+0.340 
z- - 3 ' -0.240 

1 245+ - 113 
1 - Z *' J -0.092 

1 fi52 +0 - 089 


4 lfj+0.07 
4 - iD -0.07 

4 lfi +0 - 07 
c- OO+0.28 
O-oo_0.20 

5 oq+0.28 
O.O»_ 21 

4 4n+ - 06 

4 - 4u -0.05 
4 QQ+0.06 


3.92lg;i 7 

O QO + 0.17 

°- 3o -0.27 

fi 21 + ' 56 
°- zl -0.39 

D-Jl-0.42 
4 40+0.12 
^•* a -0.09 
4 K 4 +0.04 
^•0^-0.17 


5 SQ+0-06 
°' oa -0.06 

fi 12+ 007 
D - lz -0.05 
P- 7c -+0.04 
°- '°-0.04 

C- oy + 0.05 

P- OO+0.04 
O-JJ-0.04 
c c-7+0.04 
°- ' -0.04 


14.59t°;« 
I3.99i ;i 7 
14.93l8;l! 

i4.6ii° ;il 

15.9818:13 

15.37+-0 


9 oq+0.20 
Z -J»_0.16 

1 oq+0.08 
l-OO_0.08 
9 78 +0.04 
z -'«_0.04 

2 4 5 +0.05 
z - 4o -0.09 
9 ao +0.09 
z -° z -0.09 
2 i q+0.06 
z - i9 -0.05 



Note. — Table [4] is available in its entirety in the electronic edition of the Journal. A short extract from it is shown here, for guidance regarding its form and content. Column 
descriptions: rt: model tidal radius (p(rt) = 0); i?, c : model projected core radius, at which intensity is half the central value; R^: model projected half-light, or effective, radius 
(contains half the total luminosity in projection); Rh/R c : measure of cluster concentration; Io: model central luminosity surface density in the V band; jo: logarithmic central 
luminosity volume density in the V band; Ly. total integrated model luminosity in the V band; Vtot = 4.83 — 2.5 log (Ly /Lq) + 5 log (D/10 pc): total, extinction- corrected apparent 
V-band magnitude; 1^ = Ly /2irR^ l : V-band luminosity surface density averaged over the half-light radius. Uncertainties are 68% confidence intervals, computed as described in the 
text. 

00 
03 



Table 5. Derived Dynamical Parameters for M31 Young Clusters 



Name 


Filter 


TP° P Model log M tot 


log E b 


log E 


log po 


log S h 


log <7 Pl0 


log v cscfi 


log t [h 


log fo 






[M Q L~] v ] [M ] 


[erg] 


[M@ pc" 2 ] 


[M e pc- 3 ] 


[Mq pc- 2 ] 


[km s" 1 ] 


[km s" 1 ] 


[yr] 


[M Q (pc km s" 1 )" 3 ] 



B015D F450W 0.088 



+0.01 



B015D 



B040 



F814W 0.088±o.oi 



F450W 0.094 



+0.01 



K66 
W 

K66 
W 

K66 
W 



4.83 
5.07 
4.69 
4.82 
4.30 
4.54 



+0.08 
-0.08 
+ 0.08 
0.08 
+0.06 
0.07 
+0.07 
0.07 
+0.06 
0.06 
0.06 
0.06 



48.82 
48.91 



+0.09 
-0.09 
+4.13 



-3.46 
48.76t°;°g 
48.80t°-?o 
48.25t£gg 
44 sq+ 7 - 58 



3 10 +0.09 
o,lu -0.09 
,+0.08 
-0.09 

-,+0.28 



3.1V 



4 28^ 
^• z °-0.20 

4 on + 0.29 

s+0.07 
-0.07 
,+0.07 
-0.07 



3.38] 
3.36^ 



2.86^ 

s+0.18 
-0.27 
r icr+0.56 
°' lo -0.39 
5+0.57 
-0.43 
3 46+ 12 
, +0.06 
-0.18 



2.88"* 



5.25"* 



3.51 1 



1.34 
0.78 
1.72 
1.40 
1.80 
1.16 



+0.21 
0.16 
+0.09 
0.09 
+ 0.06 
0.07 
+0.07 
0.10 
+0.10 
0.11 
+0.08 
0.07 



256+ - 039 
U - ZOD -0.042 

251 + - 076 
u - zol -0 .043 

U.28b_ 033 

290+ 031 
u.z»u_ 033 

229+ ' 031 
u.zzy_ Q 034 

n iqs+ 087 



0.914 



+0.032 
0.034 
+ 0.r,48 



924"* 
u,az -0.034 
-, n1 7 +0.036 
l.ui I _ .036 

l.U2»_o 036 
n 875 + 0.030 
U.O(O_ 32 

q 40 + - 530 



1+0.17 
-0.20 
10 61+° 12 
o 4 7 +0.04 
-0.03 

-0.06 
8 Q4+0-08 
°- M -0.08 
-,+0.05 
-0.06 



9.9V 



9.70^ 



+0.251 

-0.242 

Q15+ - 250 
u.aio_ 351 

3 096+ - 561 

°- u3u -(] .395 

i p.4+0.578 
0.101_ 437 

1 wn+ ai48 

i.o(u_ 109 
^+0.050 
0.183 



0.891 



1.604 i 



Note. — Table [5] is available in its entirety in the electronic edition of the Journal. A short extract from it is shown here, for guidance regarding its form and content. Column 
descriptions: Ty° p Ly: assumed mass-to-light ratio in the V band; Mtot = T' op Ly: integrated model mass; E b = —(1/2) Jo" 47rr 2 p<f> dr: integrated binding energy; So: central 
surface mass density; pg: central volume density; S^: surface mass density averaged over the half-light radius; <r p .o: predicted line-of-sight velocity dispersion at cluster center; v CSCi o: 
predicted central "escape" velocity; log t T ^: two-body relaxation time at model projected half-mass radius; log/p = log [po/(2 7rcr 2 ) 3 / 2 ] : a measure of the model's central phase-space 
density or relaxation time. For fo in these units, and t rc in years, logt rc ~ 8.28 — log/o jMcLauehlin fc van dcr Marcl 2005). Uncertainties are 68% confidence intervals, computed 
as described in the text. 
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